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ABSTRACT 

This routine multiplies two fixed point, single precision, 
two's complement numbers to form a double precision, two's 
complement product. 
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i - re_ t. ui re. :.::,:■ ts 

i. 1 Memory 

IK or larger alterable memory 

1.2 Equipment 

NOVA central processor 

1. 3 External Subroutines 
Unsigned multiply ( .MPYU) 

1.4 Other 
lone 

2. OPERATING PROCEDURE 

2. 1 Calling Sequence 

JSR ,MPY 
return 

2.2 Input Format 

One fixed point, single precision operand is passed 
in AC1 # the second in AC 2. 

2. 3 Output Format 

The double precision result is returned in 
AC0 (high order) and AC1 (low order). 

2.4 Error Returns 
i^one 

2.5 State of Active Registers upon Exit 

AC0, AC1, AC3 # and Carry are destroyed. AC2 
remains unchanged. 
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2.6 Cautions to User 
None 

3 . 271£H 1 Z2L 

3. 1 Algorithms 

The signed multiply routine calls the unsigned 
multiply and uses a correction factor to adjust 
the final result without resorting to a deter- 
mination of the algebraic sign of the result. 

Let Nl and N2 be the operands. Four cases must 

be examined. 

If Ml ^ 0, and K2 ^ 0, no correction is necessary. 

If Nl > and N2 < 0, then N2 as an unsigned 
number is 2**16 - abs(K2). 

The result of an unsigned multiply is thus 
2**16*abs(Ml) - abs(M)*abs(N2) . 

The true result is negative and should be 
2**32 - abs(Nl)*abs(N2). 

This can be obtained by adding 
2**32 - 2**16*abs(&l) 

to the unsigned multiply result. But since 

2**32 - 2**16*abs(Ml) = 2**16*(2**16 - abs(Ml) ) 

this is equivalent to subtracting Nl from the 
high order result of the unsigned multiply. 

If Ml < and N2 -> , a similar analysis gives 
us a subtraction of N2 from the high order 
result as the correction. 
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If ftl < and N2 < 0, the unsigned product is: 

((2**16 - abs(Ail)) * (2**16 - abs(N2) ) rood 2**32 

= -2**16*abs(Bl) - 2**16*abs(N2) + abs(fcl) *abs(N2) 

To obtain the true result of 

abs(Nl)*abs(N2) 

we merely subtract both Nl and &2 from the high 
order of the unsigned result. 

The correction factor is determined before the 
unsigned multiply is performed and subtracted 
from the high order word of the result to give 
the correct two's complement result. 

3.2 Limitations and Accuracy 
The routine is exact. 

3. 3 Size and Timing 

The routine is 16 (octal) words in length and 
requires the unsigned multiply subroutine (.MPYU). 

Average execution time is 56.4 /J seconds in 
addition to the time for the unsigned multiply. 
Average execution time for the unsigned mul- 
tiply is 340 JU seconds. Total average execu- 
tion time is thus 396.4 JU seconds. 

3.4 References 

See the unsigned multiply write-up, 093-000015. 

3 . 5 Flow Diagram s 
Bot applicable. 
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An ASCII source tape (090-000013) of .MPY is provided 
with the standard NOVA software. This tape should be 
edited into the user software that requires the signed 
multiply routine. 

5. PROGRAM LISTING 

A listing of .MPY follows. No origin is given in the 
source, enabling the user to edit the routine anywhere 
into his program. 



I MULTIPLY 

I MULTIPLIES TWO FIXED POINT* SINGLE PRECISION* 

I TWO'S COMPLEMENT NUMBERS 



I IN PUT I 



Nl IN AC1* N2 IN AC2 



I OUTPUT! 

I 



N1*N2* HIGH ORDER IN AC0# LOW ORDER IN 
AC! 



% CALLING SEQUENCE I 

I JSR .WPY 
I RETURN 



1 UNCHANGED: 
1 DESTROYED! 



AC2 

AC0# AC1* AC3* CARRY 







I REQUIRES! 


• MPYU 


CUNSIGNEO MULTIPLY) 


00000 


054813 


• MPYs 


STA 


3*.AB03 


1 


SAVE RETURN 


0000! 


176400 




SUB 


3*3 


1. 

J 


AC3 WILL CONTAIN CORRECTION 
FACTOR 


00002 


125112 




M3VL# I*1*S£C 


I 


TEST SIGN OF Nl 


00303 


157000 




ADD 


2*3 


I 


ADD N2 TO CORRECTION 


00094 


151112 




MOVL# 2*2*S2C 


I 


TEST SIGN OF N2 


00005 


137000 




ADD 


1*3 


1 


ADD Nl TO CORRECTION 


88^06 


054014 




STA 


3*.AB10 


1 


SAVE CORRECTION FACTOR 


0S§87 


006015 




JSR 


f .AB30 


I 


CALL UNSIGNED MULTIPLY 


00013 


034014 




LPA 


3*.AB1§ 


1 


GET CORRECTION 


00011 


162400 




SUB 


3*0 


1 
I 


TRUE RESULT IN 2'S COMPLEMENT 
FORM 


00012 


002013 




JMP 


».AB03 


1 


RETURN 


00013 


000030 


•A803I 







J 


SAVE AC3 


00014 


000000 


•A810S 







I 


STORAGE FOR CORRECTION FACTOR 


00015 


000015 


•AB30t 


• MPYU 


1 


UNSIGNED MULTIPLY ADDRESS 



